Axial Current Generation from Electric Field: Chiral Electric Separation Effect 
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We study a relativistic plasma containing charged chiral fermions in an external electric field. 
We show that with the presence of both vector and axial charge densities, the electric field can 
induce an axial current along its direction and thus cause chirality separation. We call it the Chiral 
Electric Separation Effect (CESE). On very general basis, we argue that the strength of CESE is 
proportional to j-ivt-^A with nv and fiA the chemical potentials for vector charge and axial charge. 
We then explicitly calculate this CESE conductivity coefficient in thermal QED at leading-log order. 
The CESE can manifest a new gapless wave mode propagating along the electric field. Potential 
observable of CESE in heavy-ion collisions is also discussed. 
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Introduction. — It was discovered and understood a 
long time ago that an external electric field, when ap- 
plied to any conducting matter, induces a vector current, 
as described by the Ohm's law 



jv = (jE, 



(1) 



where cr is the electric conductivity of the matter and 
we use the convention that the electric current is ejy. 
In quantum field theories such as the Quantum Elec- 
trodynamics (QED) and the Quantum Chromodynam- 
ics (QCD), one has not only the vector current jy but 
also the axial current when there are charged chiral 
fermions. A very interesting question is then, in addi- 
tion to the above conducting vector current under applied 
electric field, what are the other possible current gener- 
ations in response to externally applied Maxwell electric 
and/or magnetic fields. 

Recently, the QCD axial anomaly has been found 
to induce the following two phenomena in the high- 
temperature deconfined phase of QCD, the quark-gluon 
plasma (QGP), with the presence of an external magnetic 
field: the Chiral Magnetic Effect (CME) and the Chiral 
Separation Effect (CSE). The CME is the generation of 
vector current and thus the electric charge separation 
along the axis of the applied magnetic field in the pres- 
ence of nonzero axial charge density arising from fluctuat- 
ing topological charge [1-5] . With an imbalance between 
the densities of left- and right-handed quarks, parame- 
terized by an axial chemical potential /i^, an external 
magnetic field induces the vector current jy = (ipYip): 



(2) 



with chiral conductivity 0-5 = -^f. A "complemen- 
tary" effect also arising from the axial anomaly is the 
CSE which predicts the generation of an axial current, 
Ja — (^'7^75^') J and thus separation of axial charges 
along the external B field at nonzero vector charge den- 



sity (parameterized by its chemical potential /Jv) [6, 7]: 

jA = (J5t^vB. (3) 

There are robust evidences for both CME and CSE from 
kinetic theory, hydrodynamics, holographic QCD mod- 
els in strong coupling regime as well as in lattice QCD 
computations [8-19]. Experimentally the hot QGP is 
created in high energy heavy ion collisions at the Rela- 
tivistic Heavy Ion Collider (RHIC) and the Large Hadron 
Collider (LHC). In such collisions, domains of QGP with 
nonzero chirality {^ia 7^ 0) can arise from topological 
transitions in QCD and there are also extremely strong 
transient E and B fields [20-24], so the CME and CSE ef- 
fects can occur. There have been measurements of charge 
asymmetry fluctuations motivated by CME predictions 
from the STAR [25] and PHENIX [26] CoUaborations at 
RHIC as weU as from the ALICE [27] at LHC. The pre- 
cise meaning of these data is under investigations, see e.g. 
[28]. It has also been proposed that the combination of 
the CME and CSE leads to a collective excitation in QGP 
called Chiral Magnetic Wave (CMW) [29]. The CMW in- 
duces an electric quadruple of QGP that can be measured 
via elliptic flow splitting between negative/positive pions 
[30], with very supportive evidences from recent STAR 
measurements [31]. 

There is, however, one more possibility that has not 
been previously discussed, namely the generation of an 
axial current in the electric field. We find this to occur 
when the matter has both nonzero /zy and nonzero ^a- 



3a = XelJ.vfJ.AE, 



(4) 



which can be called a Chiral Electric Separation Effect 
(CESE). In this Letter, we will derive this relation, and 
explicitly compute the CESE conductivity Xe in thermal 
QED. With this new relation found, one can nicely com- 
bine all four effects into the following form: 



Jv 
3A 



a 



E 
B 



(5) 
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Chiral Electric Separation Effect. — To intuitively un- 
derstand how the CESE (4) arises, let us consider a 
conducting system with chiral fermions. When an elec- 
tric field is applied, the positively/negatively charged 
fermions will move parallel/anti-parallel to the E direc- 
tion and both contribute to the total vector current as 
in Eq.(l). If jiy > then there will be more positive 
fermions (who are moving along E) and if further /u^ > 
then there will be more right-handed fermions than the 
left-handed ones, there will thus be a net flux of right- 
handed (positive) fermions moving parallel to E. This 
picture is most transparent in the extreme situation when 
the system contains only right-handed fermions (i.e. in 
the limit of jiy = > 0), then it is obvious that there 
will be both a vector current and an axial current gen- 
erated in parallel to E. The same conclusion can be 
made when both /iy and fj.A are negative. In cases with 
/iy > > /iA or /iy < < fj,A one can follow the same 
line of argument to see an axial current generated in anti- 
parallel to the E direction. 

Different from the CME and CSE, the generation of 
axial current via CESE is not related to axial anomaly, 
but rather arising from the same conducting transport 
responsible for usual conduction in Eq.(l). To demon- 
strate that, let us consider the conduction of left-handed 
fermions, jl = cr(T; /i^, ijlr)E where = /J-v — P-A and 
^■R — fJ-v + fJ'A are the left-handed and right-handed 
chemical potentials, respectively. Due to symmetry rea- 
son, the conducting current of right-handed fermions 
must be Jr ~ cr(T; /i^, iil)E. We focus on the situation 
when all chemical potentials arc small comparing to the 
temperature T, so that cr(T; /i^, ^.r) « cro(T)-l-cro(2^)/^i+ 
a'^{T){p,\ + fij^) where the second term represents the 
chemical potential correction to the fermion conduction 
(here, the left-handed ones) and the last term arises from 
correction to the screenings (as will be seen in our ex- 
plicit calculation later). Note that due to charge conju- 
gation symmetry there will be no linear terms of fiL,R in 
the expansion. Now the total vector current would then 
be jv = JL + JR = <T{T;nL,t^R)E + a{T; fiR, hl)E w 
[2f7o + (cq + 2(TQ)(/iy + E while the axial current 
would then be Ja = Jr — Jl ~ 4(Tq fiyfiAE. We there- 
fore see that its origin is not from the axial anomaly but 
from the conduction in a chiral many-body environment, 
and the CESE is completely different in nature from re- 
cently discussed axial current generation via anomaly in 
certain superfluid systems [32, 33]. 

Computation of the CESE Conductivity. — We will now 
compute explicitly the leading-log order CESE conduc- 
tivity for thermal QED plasma using the Kubo formula 
under the condition fj,v,fJ-A <C T. The extension of this 
computation to QCD is straightforward, and will be pre- 
sented elsewhere. 

Let us denote — XefJ-vfJ-A- Starting with the re- 



tarded vector-vector and vector-axial correlators G 

iRij 



Rij 
VV 



and (Te are given via the Kubo formulae as: 

3 

a = V lim hm Ag^»^(^ fc), (6) 

i=l 
3 

'^^ = 5]hm hm -f G2t)(^,fc), (7) 

In Fig. 1, the shaded circle represents an effective ver- 
tex, see Fig. 2. This effective vertex represents a resum- 
mation of a set of ladder diagrams which contribute to 
the same order owing to the pinching singularity when 
a;,fe — [34]. The dotted line in Fig. 2 represents the 
hard thermal loop (HTL) resummed propagator. All 
other kinds of diagrams (e.g., the box diagrams which 
are vanishing identically due to Furry theorem at zero 
chemical potentials but finite at nonzero chemical poten- 
tials) are sub-leading-log order. 




FIG. 1: The Feynmann diagrams for retarded correlators 
G^v (Left) and G^y (Right). 

One can write down the retarded correlators explicitly: 
G'^y = -eJ^TT[r\P + K,P)SiPh^S{P + K)], (8) 
G'^A = ~e j^Tr[r{P + K,P)S{PW^,S{P + K)\{%) 
The effective vortex in the above is given by Fig. 2, 
T^'{P + K, P)=Y + e^ [ YS{P + K + Q) 

Jq 

xr^(P + K + Q,P + Q)S{P + Q)-l^*Dp,{Q). (10) 

The *Dpa{Q) is the HTL propagator for photon while 
S{P) is the electron propagator at nonzero /iy and /iA, 

S{P) 



^ -7 • (Ps - Ss) 



s=± 



{Ps - ^sY 



(11) 



where S± = (1 ± 75)/2 is the chirality projection, P!^ = 
{P^,p) with P^ = Po + M± and /i± = /iy ± ha. Sub- 
stituting all these into the Kubo formulae and following 
essentially the steps as in Refs. [35, 36], we finally obtain. 



s,a— ± 

CTe = ^ Sasa, 
s,a— lb 



(12) 
(13) 



and (as diagrammatically shown in Fig. 1), the a 



-a- I ■7:r^n'p{p- aiJ.s)xl{p), (14) 



3 J (27r)3 
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with np{x) = l/[exp(a;/T) + 1] and xiip) = 
'^'aip)/K- The is defined through 2p^V%{p) ^ 
Us{p)TX{p)u,{p), 2pWt{p) = V-si-p)rUp)v-si-p), 
with Fj. = — 2ImI]^g(pg = ±Ep,p) the decay width for 
fcrmions (+) and anti-fermions (— ) of chirahty s. 





FIG. 2: Integral equation for the effective vertex F**. 

The vertex integral equation (10) can be finally recast 
into a differential equation for Xaip)- 



2am'^ ln(l/e) 



+ 



aTmj^ ln(l/e) 



P^ 
2T^ 



X'aip) 



2nF[p- a^s])) P 



p^ <Px'L(p) 

2 dp'^ 

d xlip) 
dp 



(15) 



where nsix) = l/[exp(x/r) — 1] is the Bose-Einstein 
function, a = e^/in, m'j. ~ ^ (t^ + is the effec- 
tive mass of electrons, m|, — ^ X]s=± + 
the Debye screening mass of photons. We note that the 
leading-order terms in the effective vertex and those in 
the fermion decay width cancel each other and as a result 
the leading-log results for conductivities are actually at 
the order Cl(l/[e3 ln(l/e)]) rather than 6(l/[e ln(l/e)]) 
as expected from asa ~ e/F^. 

Following Ref. [37], we solve this differential equation 
numerically by using variational method. This varia- 
tional scheme converges very fast, and at very high pre- 
cision we obtain (up to 0{^s/TY order) 



T 



e3 ln(l/e) 



-0.62224 



3.9238- 



rp2 



4.59867a^ -I- 2.56237^ 



(16) 



and thus 



T 



e3 ln(l/e) 
20.499-^^^-^ 



15.6952-^ 7.76052 
T 



2^2 



T2 e3ln(l/e)' 



,(17) 
(18) 



aforementioned CMW [29, 30] in which the fluctuations 
of vector and axial currents are coupled together by ex- 
ternal B field to form a propagating wave. Now the new 
CESE effect introduces nonlinearity (aka the fJ-vfJ-A term) 
and makes the problem more nontrivial. 

Let us consider a thermal QED or QCD plasma in the 
static and homogeneous external E, B fields, and study 
the coupled evolution of the small fluctuations in vector 
and axial charge densities. The presence of vector den- 
sity and current will induce additional electromagnetic 
fields so that Etot ^ E + 6E and Btot ^ B + 5B with 
6E (X ejy and 5B cx ejv- As in the case of CMW, one 
can first replace the chemical potentials in Eq.(5) with 
the corresponding charge densities: — ctv,Ajv a^ 

where the ay, a are the susceptibilities defined as ay, a = 
diiv,A/djy j!^. These relations are valid as long as the 
chemical potentials are small compared with temperature 
T. Then, combining Eq.(5) with the currents' continuity 
equations dtjyj^ + V ' Jv.A = and Maxwell's equation 
V • Etot — ejy and V ' -Btot = 0, one can obtain: 

dtjy + eaoiv + <^5aA{B ■ v)^^ 
+2a2a^vJviE ■ v)jv + '2'^2al f^E ■ v)Ja = , 
ddA + <^5av{B ■ v)i° + XeavaAjviE ■ v)j° 
+XeayaAjl(-E- V)i° =0. (19) 

where we have introduced ctq and CT2 defined through the 
small chemical potential expansion of ct, cr = co+o'2(A'y + 
/i^). In arriving at the above, we have made a few ap- 
proximations: first we keep only terms up to O(jy^) or- 
der (such that terms like j^SE-xjj^ ^ 0{j^) and (y2l^'^jv 
are dropped); second we assume external fields are ex- 
tremely strong such that the magnetic field feedback 
terms ^ a^a6B - ^z ^ a^aeiyjdj) are negligible to other 
0{p) terms ^ Xeci^E{jdj) and ^ a2a'^E{jdj). We wiU 
next linearize the above equations by considering fluctu- 
ations on top of small uniform background vector and 
axial densities ny and ua, so that jy^ — nv,A + ^Jva 
with Sjyj^ <^ ny.A being fluctuations. Note that nv,A 
themselves must still be small (compared to T^) to en- 
sure the linear relations between jy^ and ^v,A- Strictly 
speaking only a uniform axial density ha is a static solu- 
tion of the above equations while a uniform vector density 
ny suffers from the damping term eaojy and is only ap- 
proximately static on time scale shorter compared with 



1/euQ. Keeping only linear terms in Sjy j\ we obtain: 



In the limit — > 0, our result for a is in agreement with 
(7 w 15.6964 -i 1 T, ; ^ obtained in Ref. [37]. 

ln(l/e) L J 

Coupled Evolution of the Two Currents. — As seen 
in Eq.(5), with the presence of external electromagnetic 
fields, the vector and axial currents mutually induce each 
other and get entangled in a nontrivial way. It is of great 
interest to understand the coupled evolution of small fluc- 
tuations of the two currents. A very good example is the 



dtSjl 4- e<TQ5jl + a5aA{B ■ v)<5jA 

+2<T2a'^ny{E ■ \7)6j^ + 2a2a\nA{E ■ \7)SjA = , 
dtSjA + o^aviB ■ \7)5jy + Xe^vaAnviE ■ vJ^Ja 



+XeavaAnA{E ■ v)^i° = • 



(20) 



In order to find possible normal modes, we make the 
Fourier transformation of these equations. Using Sjy = 



4 



L k k) e-*('^*-'= =«), we obtain from Eqs. (20) the 

following relations: 

uCv + ieaoCv — a^aAiB ■ k)CA 

-2a2alrnv{E ■ k)Cv - 2a2a\nA{E ■ k)CA = 0, 
(jjCa — <750iv{B ■ k)Cv — XeCtvOiAnA{E ■ k)Cv 

-XeayoiAnviE ■ k)CA ^ ■ (21) 

Without loss of generality, we can always assume B is 
along z-axis, i.e., B = Bz while E — Ee. The dispersion 
relation obtained from Eq. (21) can be expressed as 

u; = [zecTQ - v+{e ■ k)] 



±^\J[ieao~v.{e-k)Y +^A^{k), (22) 

where v± = ± Va with = 2<j2aynv E and Va — 
XeCtvO-AnyE, and 

Ax{k) = [(75aAB{z ■ k) + 2a2aA'>T'AE{e ■ k)] 

X [a5avB{z ■ k) + XeCtvCiAnAE{e ■ k)] .(23) 

To manifest the physical meaning of the solutions in (22), 
let us consider two special cases. 

(1) The case with only B = Bz and E ^ 0. 

Then Eq.(22) reduces to w = ±y/{vxk^)^ - {eao/2)'^ - 
i{eao/2) with speed — <J5y/avciAB. When Vy,kz ^ 
eao/2 we get two well-defined propagating modes w « 
ztv^kz — i{eao/2): these are generalized CMWs which 
reduce to the CMW in [29] when ao — and ay = 
aA- When v^k^ < ef7o/2 the two modes become purely 
damped. 

(2) The case with only E ^ Ez and B = 0. 

First, we consider a background without vector density, 
i.e., nv — 0. In this case, we find two modes from (22) 



LO = ±^{vekz)^ - (eao/2)2 - i(eao/2) 



(24) 



with Ve = aAnA\J2a2XeOLvOLAE. Similar to the CMWs, 
when Vekz 3> ecro/2 there are two well-defined modes 
bj w iwefc^ — i(e(To/2) from CESE that propagate along E 
field and can be called the Chiral Electric Waves (CEWs). 
They become damped when vji^ ^ eCTo/2. 
Second, if the background contains no axial density, i.e., 
ua = 0, then we see that the vector and axial modes 
become decoupled, and Eq. (22) leads to 



^v(k) 

UJA{k) 



Vykz - iieao) 

Vakz- 



(25) 



The first solution ajy(fc) represents a "vector density 
wave" (VDW) with speed Vy — 2(j2oiynvE that trans- 
ports vector charges along E field but will be damped 
on a time scale ~ l/(eCTo). The second solution 0JA[k) is 
a new mode arising from CESE and represents a propa- 
gating "axial density wave" (ADW) along E with speed 
Va = XeOivoiAnvE and without damping. 



Summary and Discussions. — In summary, we have 
found a new mechanism for the generation of axial cur- 
rent by external electric field in a conducting matter with 
nonzero vector and axial charge densities, which we call 
the Chiral Electric Separation Effect. We have computed 
the CESE conductivity coefficient in a QED plasma and 
also studied possible collective mode arising from it. 
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FIG. 3: A schematic illustration for CESE induced net charge 
distribution and correlation patterns in Cu + Au collision. 

We end by discussing possible observable effect that 
may be induced by CESE in heavy ion collisions. In 
the created hot QGP there can be both vector and axial 
charge densities from fluctuations and topological transi- 
tions. There are also very strong electric fields during the 
early moments of heavy ion collisions [20-22, 24]. One 
particularly interesting situation is in the Cu + Au col- 
lisions, where due to the asymmetric nuclei (rather than 
from fluctuations) there will be a strong E field direct- 
ing from Au nucleus to Cu nucleus [40] . In this case the 
E field will lead to both an in-plane charge separation 
via (1) and an in-plane chirality separation via (4). The 
resulting in-plane axial dipole will then further separate 
charges via CME along the magnetic field in the out-of- 
plane direction, and cause an approximate quadrupole 
at certain angle q in between in- and out-of-plane . We 
therefore expect a highly nontrivial charge azimuthal dis- 
tribution pattern SN±{(j)) ^ dcos{(j)) + qcos{2(j) — 2'^g) 
with the dipole term due to usual conductivity and the 
quadrupole term due to CESE and CME effects. This 
pattern may possibly be measured either via charged pair 
correlations or the charged multiple analysis [22]. More 
quantitative predictions will require proper modeling of 
the QGP and solving the Eq. (20) which will be reported 
in a future work. 
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